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ABSTRACT 

We describe an upgrade to the Cosmic Background Imager instrument to increase its sur- 
face brightness sensitivity at small angular scales. The upgrade consisted of replacing the 
thirteen 0.9-m antennas with 1.4-m antennas incorporating a novel combination of design 
features, which provided excellent sidelobe and spillover performance for low manufacturing 
cost. Off-the-shelf spun primaries were used, and the secondary mirrors were oversized and 
shaped relative to a standard Cassegrain in order to provide an optimum compromise between 
aperture efficiency and low spillover lobes. Low-order distortions in the primary mirrors were 
compensated for by custom machining of the secondary mirrors. The secondaries were sup- 
ported on a transparent dielectric foam cone to minimize scattering. The antennas were tested 
in the complete instrument, and the beam shape and spillover noise contributions were as 
expected. We demonstrate the performance of the telescope and the inter-calibration with 
the previous system using observations of the Sunyaev-Zel'dovich effect in the cluster Abell 
1689. The enhanced instrument has been used to study the cosmic microwave background, 
the Sunyaev-Zel'dovich effect and diffuse Galactic emission. 

Key words: instrumentation: interferometers - methods: data analysis - cosmic microwave 
background - X-rays: galaxies: clusters. 



1 INTRODUCTION 

The Cosmic Background Imager (CBI) JPadin et alj|2002l) was a 
13-element co-mounted interferometer operating at 26-36 GHz, 
designed primarily to observe the power spectrum of fluctuations 
in the cosmic microwave background (CMB) on angular scales 
of 5 arcmin to 1 deg (multipoles I ~ 400 to £ ~ 3500). Between 
2000 lanuary and 2005 April, the CBI operated from the Cha- 
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jnantor Plateau, Chile at an altitude of 5100 m and during this 
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200 1 l2009h . In addition it was also used to make observations 



of the Sunyaev-Z erdovich (SZ) effect in a sample of low-redshift 
(z ^ 0.1) clusters dUdomprasert et all] 20041) . and measurements of 
'anomalous' microwave emission from dust in a range of Galac- 
tic objects |Casassus et ail |2004 120061. 120081 ; lHales et alj |2004| : 
iDickinson et al.ll2006l . l2007h . 
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Figure 1. Left: The change in thermodynamic temperature relative to the CMB due to the thermal SZ effect from three clusters, at z = 0.15, with different 
large-scale properties. The vertical lines represent radii at r2soo (dashed-dotted), r$oQ (dashed) and rjoo (solid). The profiles are (red: isothermal beta model) 
A7\ = A7"o(l + (r/n) 2 )'/ 2-3 ' 3 / 2 with r\ = 0.75arcmin, j3 = 0.7 and A7o = l.OmK; (blue: suppressed large scales) A72 = A7) exp(— {r/riY) with ri = 7 arcmin; 
and (green: additional large scale component) A73 = A7i + 0. 2exp(— (r/ri) 1 ) mK. Right: The corresponding visibility amplitude with 1 .4-m diameter antennas 
as a function of baseline length u, at V = 3 1 GHz. The vertical lines represent baseline lengths corresponding to the three radii (e.g., 112500 = 1 / '"250o)- The error 
bars represent typical data from the CBI2 experiment, including both thermal noise and intrinsic CMB components. 



These observations were made using antennas 90 cm in diam- 
eter. In 2005 - 2006, the CBI was upgraded to larger 1 ,4-m antennas 
('CBI2') to increase the effective collecting area and to allow obser- 
vations at higher resolution without compromising surface bright- 
ness sensitivity. Observations with the CBI2 continued until 2008 
June, afte r which its site and mount were u sed for the QUIET ex- 
periment dOUIET Collaboration et alj201Ch . During this period the 
CBI2 completed a programme of observations of diffuse Galac- 
tic emission, the CMB powe r spectrum and targeted SZ clusters 
dPickinson et alj I2009L I20T0I ; ICastellanos et alj|201 ll : Ividal et~afl 
|201ll and further papers in preparation). In this paper we describe 
the antenna design that was used in the CBI2 upgrade. We sum- 
marize the main science goals for the upgrade and present com- 
missioning results that confirm its effectiveness. We also present a 
combined analysis of an SZ detection in the cluster A 1689. This 
cluster was observed both with the original CBI (hereafter 'CBI1') 
and with the upgraded CBI2 and allows us to demonstrate both the 
inter-calibration of the two instruments and the benefit of measur- 
ing the SZ decrement with the larger CBI2 antennas. 



2 SCIENCE MOTIVATION 

The angular scales to which an interferometer is sensitive are set 
by the lengths of the baselines between the antennas, with longer 
baselines responding to finer-scale information in the sky bright- 
ness. However, for a fixed antenna size, the sensitivity of a baseline 
to extended sources decreases rapidly as the baseline is lengthened. 
In the Rayleigh- Jeans limit, the temperature sensitivity AT is given 
approximately by 

AT = X 2 AS/{2kfQ), (1) 

where A5 is the flux density (point source) sensitivity, is the solid 
angle of the main lobe of the synthesized beam, and the filling fac- 
tor, /, is the fraction of the synthesized aperture that is filled with 
antennas. This is simply a modification of the Rayleigh- Jeans equa- 
tion to reflect the fraction of photons captured instantaneously by 
the aperture - the exact temperature sensitivity as a function of an- 
gular scale will depend of the configuration of the antennas within 
the synthesized aperture. Increasing the resolution of an interfer- 



ometer without losing brightness sensitivity thus requires that ei- 
ther the number of antennas be increased, or the antenna size be 
increased, in order to maintain the filled fraction of the synthesized 
aperture. If the number of baselines is fixed, and the antennas are 
not changed, lengthening the baselines results in an increase in inte- 
gration time to reach the same temperature sensitivity proportional 
to the fourth power of the baseline length. 

The primary goal of the CBI2 upgrade was to increase the 
temperature sensitivity of the instrument on its longer baselines of 
3 - 5.5 m, i.e., corresponding to angular scales of 6 - 12 arcmin, 
on which the CBI1 array was not well filled. Improved sensitivity 
on these longer baselines would provide significantly improved ob- 
servations of the SZ effect in massive galaxy clusters. In CBI1 SZ 
observations, the shortest baselines were heavily contaminated by 
primary CMB anisotropies, while the longer baselines lacked ther- 
mal sensitivity. Moderately massive clusters typically have virial 
radii of ~2Mpc, which at a redshift of z ~ 0.15 corresponds to 
an angular size of ~ 12 arcmin. This is well-matched to the new 
CBI2 array, which is thus able to measure the cluster gas out to the 
outskirts of the clusters with significantly less contamination from 
primary CMB fluctuations than was the case for CBI1. 

The motivation to concentrate on measuring the SZ effect out 
to the virial radius in complete samples of clusters was driven by 
the need to further understand the X-ray-SZ and weak-lensing-SZ 
scaling relations in support of SZ survey experiments. The SZ ef- 
fect measures the Comptonization parameter, y = J kT e / (m e c 2 )dZ, 
which is proportional to the electron pressure integrated along the 
line of sight. SZ surveys are designed to measure the integrated SZ 
effect, Y = JydSl, providing empirical measurements of the clus- 
ter co-moving SZ luminosity function dN/dY. However, in order 
to relate these measurements to cosmology via the cluster mass 
function, AN /AM, a well-calibrated relationship between Y and 
the total mass M is required. This can be achieved by combining 
SZ measurements of known clusters with X-ray and weak lens- 
ing data, along with modelling that accurately describes the dis- 
tribution of the cluster components (dark matter, gas and galax- 
ies) in a way that can be constrained by the observational data. 
There have been a number of recent measurements of the scaling 
between t he integrated SZ effect and the total mass, f rom both hy- 
drostatic jBenson et alJl2004l ; lBonamente et alj|2008h and gravita- 
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tional lensing jMarroneetail2 009) estimates. However these rela- 
tionships have generally only been obtained out to relatively small 
radii (~ 200 — 400 k pc) and observations with experim ents such 
as CB I2, APEX-SZ dSchwan et al1l2003l) and AMIBA JHq et alj 
1 20091) are expected to provide constraints out to a few Mpc (~ r200> 
where r x denotes the radius within which the average density is x 
times the critical density). 

FigureQ] illustrates why measurements at large angular scales 
relative to the core of the cluster are important in determining 
true SZ profiles. It shows the thermal SZ effect for toy models of 
three clusters at z = 0.15, with similar cores but different large- 
scale properties. The red curve is a standa rd isothermal beta model 
dCavaliere & Fusco-Femiaiiolll976l.[T973) . 

Ar = Ar (l + (r/r core ) 2 ) 1 / 2 - 3 / i / 2 (2) 

where ATq is the central temperature decrement, r cole is the cluster 
angular core radius and j8 controls the shape of the profile. The 
blue profile has the same radial behaviour within the cluster centre 
but has an exponential decline which becomes significant beyond 
r > r25oo (corresponding, for example, to a declining temperature 
profile). The green profile has a similar radial behaviour to the blue, 
but with an additional additive large-scale component that produces 
a larger overall value for the central SZ signal. The corresponding 
visibility profiles (as a function of baseline length in wavelengths 
w) are obtained by taking the amplitude of the fourier transform of 
the product of the SZ model, AT(r), with the primary beam of the 
interferometer B(r) (here assumed to be a Gaussian of half-power 
width 30 arcmin). The signal is converted from temperature to flux 
density units using the Planck equation. For a circularly symmetric 
model this is easiest to implement using a Hankel transform, 

2k x 2 e x f°° 
A5(«) = - ? ._ r -~ 2w AT(r)B(r)J (27Cru)rdr, (3) 

c z (e x — \y ./() 

where x = hv/kTcMB an d Jo is the zeroth-order Bessel function. 

These visibility profiles show that interferometric experiments 
cannot distinguish between the SZ profiles on baselines greater 
than 400 A. (The same would be true for total power measurements 
where the data are spatially filtered on scales greater than the equiv- 
alent angular scale, here about 10 arcmin.) However the integrated 
SZ flux density, which corresponds to the total thermal energy in 
the cluster and is the quantity which is measured as a proxy for 
mass in cluster surveys, varies by almost a factor of two between 
these cases. It is therefore important to observe on baselines short 
enough to distinguish between different large-scale cluster proper- 
ties. In the case of observing frequencies ~ 30 GHz, this requires 
an interferometer with baselines smaller than around 4 m. 

FigureQ] also displays the expected data that would be ob- 
tained from CBI2 observations of these three SZ profiles. The er- 
ror bars include both a thermal noise component and a component 
due to the intrinsic CMB anisotropy, the latter being significant on 
baselines shorter than 250A. The uncertainty due to the thermal 
noise integrates down with the square root of the observing time, 
whereas the contamination due to the primordial CMB fluctuations 
does not. Without additional frequency information which can dis- 
tinguish between the CMB and the characteristic spectrum of the 
SZ effect, sensitivity to the large-scale SZ effect is ultimately lim- 
ited by the primordial CMB fluctuations. The CBI2 baselines how- 
ever provide an excellent compromise between primary CMB con- 
tamination on the one hand and resolving out of the largest scale 
emission on the other. 



3 UPGRADE OF THE CBI ANTENNAS 
3.1 CBI1 antenna design 

The original antenna design for the CBI1 was an on-axis Cassegrain 
with a 0.9-m diameter, / = 0.33 primary, and a 155 m m-diameter 
hyper boloidal secondary with eccentricity of 1.41 ( Padin et alj 
20021) . The secondary was supported on a transparent polystyrene 
quadrupod, and the whole antenna enclosed in a can rising to 
400 mm above the rim of the primary. This can was designed to 
reduce the coupling from the secondary of one antenna to the feed 
of the adjacent antenna, and was measured to reduce such coupling 
from —90 dB to ~ —120 dB. The secondary mirror was over- 
sized, in the sense that it extended beyond the radius required to 
reflect a ray from the feed to the edge of the primary. This is a 
common feature of Cassegrain designs, and is intended to increase 
the aperture efficiency. By increasing the size of the secondary, the 
diffraction beam of the secondary is reduced, resulting in a side- 
lobe which would otherwise have missed the primary edge hitting 
the primary, and therefore contributing to the main aperture illumi- 
nation. It also has the effect however, of providing additional direct 
ray paths from the secondary over the edge of the primary, in the 
direction of the adjacent antennas. To alleviate this problem, the 
original antennas were provided with the shield can, which greatly 
reduced the spillover in the backward direction, and redirected the 
spillover power to the sky in the general direction of the main beam. 



3.2 CBI2 antenna design 

The new antenna design for CBI2 was intended to make maximum 
use of the physical area of the platform on which the CBI1 antennas 
were co-mounted. The largest antenna size that could be accommo- 
dated on the table while still using all 13 antennas was 1.4 m diam- 
eter. The design had to be cost-effective and reasonably quick to 
implement, and so was based on a commercially available reflector 
with a nominal diameter of 1.37 m and actual maximum diameter 
1.41 m (the difference being due to the roll-off of the surface at 
the rim). The focal length was 457 mm, giving an /-ratio of 0.33, 
very similar to the original CBI1 design. The reflector was fabri- 
cated by spinning, in which a circular aluminium sheet is pressed 
over a spinning mould with a roller. The surface is then rolled over 
a circular tube at the rim, and finally a circular tube is riveted to 
the back of the dish at 300 mm radius to provide a mounting point. 
This method is very quick and cheap, and measurements with a co- 
ordinate measuring machine showed that the surface accuracy 8x 
of a sample dish was better than 0.2 mm rms on small scales over 
most of the dish surface. However, the dishes also typically had 
large-scale distortion of the form 8x °= r 2 cos(20), consistent with 
the rim tube being elliptical with a deviation from circular of sev- 
eral millimetres, forcing the surface out of its paraboloidal shape. 
This error was dealt with by modifying the secondary optics as de- 
scribed below. 



3.3 Optical design of the CBI2 secondary optics 

The main competing design drivers were aperture efficiency versus 
sidelobe spillover. The existing CBI1 feed horn was modelled using 
the CorrucQ software package and the resulting feed pattern used 
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Figure 2. Ray diagrams of the CBI2 antenna with an oversized secondary 
with (top) no reshaping, showing the potential for spillover past the pri- 
mary and (bottom) with the secondary reshaped at large radii to redirect the 
spillover rays back on to the primary. 

to illuminate a model of the primary dish in the GRASPS0 soft- 
ware package. GRASP9 enables full physical optics plus physical 
theory of diffraction simulations to be done on the complete op- 
tical system. This method takes into account the fields from both 
the surface and edge currents on the reflectors, as well as blockage 
and multiple reflections (e.g., in the region of the primary shad- 
owed by the secondary). In order to minimize spillover without 
sacrificing too much aperture efficiency, the CBI2 optics design in- 
corpora ted a secondary m irror which was both oversized and re- 
shaped l lHoller et alj|2008h . In order to find the optimum balance 
between aperture efficiency and sidelobe spillover, the secondary 
mirror size was increased from the nominal ray optics size, in- 
creasing the aperture efficiency due to diffraction effects, until the 
increasing blockage began to reduce the efficiency again. The edge 
of the secondary was then reshaped by adding a quadratic term to 
the hyperboloid, starting at a point near the ray optics illumination 
edge (i.e., the point where an on-axis ray striking the edge of the 
primary would strike the secondary). This has the effect of direct- 
ing radiation closer inwards on the primary than would otherwise 
be the case. The ray traces are shown in Figure [2] (ray traces are 
not accurate modelling tools for cases such as this where the an- 
tenna properties are dominated by diffraction effects, but are useful 
to visualize the design concepts). Both the starting point and am- 
plitude of the quadratic term were adjusted to achieve a reasonable 
compromise between aperture efficiency and sidelobe level. In ad- 
dition, the distance between the primary and secondary mirrors was 
adjusted from the geometric optics value in order to maximize the 
forward gain, to take account of the fact that the optics is all in the 
near field of the feedhorn. This resulted in a shift of the secondary 
position by 5 mm towards the primary and an improvement of the 
forward gain by about 25 per cent. 

The central 14 mm of the secondary was reshaped into a cone 
of semi-angle 86 deg, such that rays from the feed close to the axis 
are not reflected directly back on to the cryostat window, causing 
standing waves. Figure|3]shows cuts in one plane through the calcu- 



2 TICRA: http://www.ticra.com/what-we-do/software-descriptions/gn 
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Figure 3. Comparison of calculated beam patterns of the CBI1 and CBI2 
antennas across the observing band of 26-36 GHz. Top: Far-out beam pat- 
terns: Black, 1.4-m antenna at 26 GHz; blue, 1.4-m at 31 GHz; green 1.4-m 
at 36 GHz; red 0.9-m antenna at 3 1 GHz (no shield can). Note that the main 
spillover lobe at 100 deg is lower in all cases for the larger antenna and 
is negligible at the top of the observing band. Bottom: Main lobe of the 
calculated beam patterns. Solid lines indicate the 1.4-m antenna. Dashed 
lines indicate the 0.9-m antenna. In both cases the colour scheme is green, 
36 GHz; red, 31 GHz and black, 26 GHz. 

lated beam patterns at 26, 31 and 36 GHz. Also shown for compar- 
ison is the calculated pattern for the 0.9-m antennas (without shield 
can) at 31 GHz. It can be seen that the spillover lobes are reduced 
by around 20 dB compared to the previous design, and that there is 
essentially no spillover lobe at the top end of the observing band. 

The main beam plots in Figure [3] also show that the forward 
gain of the CBI2 antenna is 3 dB greater at each frequency com- 
pared to the old CBI1 design, i.e., the effective aperture is bigger 
by a factor of two. This is smaller than the nominal area increase 
((1.4/0.9) 2 = 2.4) due to the under-illumination (or steeper edge 
taper) caused by the reshaping, which was necessary in order to 
keep the spillover lobes to a minimum. 

3.4 Correction for primary asymmetry 

The optical design was developed using an ideal model of the 1.4 m 
primary dish. However, as described earlier, the manufacturing pro- 
cess used to make the primary dishes introduces a large-scale de- 
formation, as a result of using a non-perfectly circular reinforcing 
ring at the rim of the dish. This results in a non-circular beam- 



The Cosmic Background Imager 2 5 




pattern with reduced forward gain. To compensate for this effect 
the surface profile of each CBI2 primary dish was measured along 
several circular tracks at differing radii. For all but one of the dishes 
it was possible to get a good fit to the distortion using just the 
second-order Zernicke polynomial Z? = r 2 cos(20) (the other dish 
also required an inclusion of a Zj = r 2 cos(30) term). By adding 
the appropriate Zernicke polynomial to each of the secondary an- 
tenna profiles it is possi ble to cancel out the effect of the distortion 
to'Sullivanet4]l2008h . This was done for each dish in turn, us- 
ing the ZEMA>U optics modelling package in order to determine 
the amplitude of the polynomial correction needed to maximize 
the Strehl ratio in a ray optics simulation. GRASP9 simulations 
of each of the antennas in turn were made using the appropriately 
corrected secondary to verify the design. The tolerance of the beam 
shape to focus position was also modelled, as the corrected optics 
were noticeably less tolerant to focussing errors than the ideal op- 
tics. The resulting secondary mirror designs (oversized, reshaped 
and with an appropriate Zernike polynomial added) were then ma- 
chined from solid aluminium using a CNC milling machine, with 
care taken to indicate the axis of the polynomial on the secondary 
such that it could be aligned with that of the corresponding primary 
dish. 



3.5 Antenna assembly 

The new antennas were mated to the existing CBI1 receivers using 
an existing mounting plate that in the original design mounted the 
receiver, primary mirror and shield can. As in the original design, 
in order to avoid introducing scattering in the beam, the secondary 
mirror was supported using a transparent dielectric material. We 
used a hollow cone of Plastezot^]LD45, an expanded low-density 
polyethylene foam (the '45' refers to the density in kgm~ 3 ). Sam- 
ples of the LD45 were tested in the lab for both their thermal and 
electrical properties. Dielectric loss was unmeasurably small at 30 
GHz - scaling from the volume fraction of the foam the expected 
value of the loss tangent would be tan 8 = 2 x 1CP 5 - and the di- 
electric constant was measured as e r = 1.06, which is equal to the 
dielectric constant of solid polyethylene diluted by the solid vol- 
ume fraction of the foam. The coefficient of thermal expansion was 



3 Zemax Development Corporation, http://www.zemax. com/| 

4 Zotefoams pic, http://www.zotefoams.com/pages/en/datasheets/ld45.htm 



however significant at about 5 x 10~ 5 K . Given a possible tem- 
perature range of > 20 K on site and the height of the cone of 400 
mm, the resulting change in focus position approaches the maxi- 
mum tolerance of the design, at around ±0.5 mm. Care was there- 
fore taken on assembly to ensure that the nominal focus position 
was achieved at the mid-range of expected ambient temperatures 
(about —5° C). Each secondary mirror was attached to a foam lid 
using a metal plate screwed in to the back of the mirror, and ma- 
chined alignment jigs were used to hold the mirror in place while 
the lid was glued to the cone. The cone was assembled from sec- 
tions bandsawed out of 100-mm thick LD45 sheet, glued together 
on joints perpendicular to the optical axis of the antenna using the 
impact adhesive Evostik TX52^f]. The adhesive joints were very 
much thinner than a wavelength, and tests in waveguide showed 
that they had negligible attenuation or reflection at 30 GHz. A dia- 
gram of the antenna assembly is shown in Figure [4] 



4 RE-COMMISSIONING TESTS 

All thirteen CBI2 antennas were assembled on-site in Chile prior 
to mounting on the CBI platform (Figure [4]l. The pointing of each 
antenna was assessed by making 5-point observations of bright cal- 
ibrators such as Jupiter and Tau A. Here the instrument is first 
pointed on-source, and then off-source at 4 pointings where the 
amplitude should be half of the total signal from the source. Since 
the CBI2 is a co-mounted interferometer the pointing errors associ- 
ated with any individual antenna have to be separated by modelling 
the response of each baseline to the combined pointing error of its 
pair of antennas and solving for the individual antenna pointing er- 
rors. These differential pointing errors were corrected for by plac- 
ing shims under the relevant mounting feet of each antenna. The 
residual individual pointing errors after this process were typically 
~ 0.5 arcmin. 

The primary beam of the new system was measured using ob- 
servations of Jupiter on a grid of 1 1 x 1 1 pointing centres spaced 
by 7 arcmin (i.e., covering offset positions of ± 35 arcmin in az- 
imuth and elevation). The integration time per pointing was 45 s. 
The resulting beam patterns for one of the CBI2 antennas both with 
and without a corrected secondary mirror are shown in Figure [5] 
and clearly show the improvement in the circularity of the beam 

5 Bostik, http ://www.bostik.co.U k/diy/product/evo-stik/TX528/9 
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Figure 5. The effect of the correction for primary astigmatism on the beam shape of the CBI2 antennas. Left: Beam pattern measured from an observation of 
Jupiter in a single, 1 GHz frequency channel centred at 31.5 GHz for a single antenna with typical primary astigmatism but fitted with a symmetric subreflector. 
Middle: The beam pattern from the same antenna fitted with its individually corrected subreflector. Contours are 5, 10, 20 ... 90 per cent of the beam peak. 
Individual antenna patterns are solved for from the visibility data for the whole array as it is scanned over the source. Right: The measured profile of the CBI2 
beam at 31.5 GHz. The scatter points are the measured data from all the CBI2 baselines in a single, 1 GHz frequency channel centred at 31.5 GHz. The solid 
line is the beam simulated in GRASP9, and the dashed line is a Gaussian of FWHM 27.8 arcmin (the same FWHM as the calculated beam), which fits the 
beam well within the half-power points. 



when using the corrected secondary. Figure [5] also shows the mea- 
sured radial profile of the beam at 31.5 GHz along with the simu- 
lated GRASP9 beam. The measured beam can be fitted to the half- 
power points with a Gaussian model with a FWHM of 27.8 arcmin 
at 31.5 GHz, also shown in Figure(5] 

To calculate the expected thermal noise of the new CBI2 array 
we assume an effective antenna collecting area of 0.8 m 2 , effec- 
tive bandwidth per channel of 0.85 GHz, correlator accumulation 
time of 4.2 s, a nominal system temperature of 30 K and a system 
efficiency (due to non-flat passbands, phase errors etc.) of 90 per 
cent. This gives an expected rms thermal noise of 1.9 Jy per sam- 
ple, or 3.9 Jy s 1 / 2 . In order to measure the actual sensitivity a series 
of blank field observations was made. To reduce the contribution 
from ground spill the observations were taken in differenced mode 
i.e., with a lead and trail field observed at the same declination as 
the main field but separated by 8 minutes in right ascension. The 
resulting sensitivity was then found by calculating the mean rms of 
the real and imaginary parts of the visibilities after subtraction of 
the lead and trail fields. Figure[6]shows a histogram of the measured 
noise from CBI2 observations of blank fields over the period 2007 
April 1 to 2008 May 10. Each datum was generated from 23 sam- 
ples of the correlator output. The measured peak in the rms noise 
histogram lies at the expected value, but with a tail in the distri- 
bution to larger noise values due to baselines containing receivers 
with higher than nominal system temperatures. 



5 SZ DETECTION OF A1689 AND COMPARISON WITH 
CBI1 

As a final check on the effectiveness of the CBI2 upgrade 
we present an analysis of an SZ de t ection in the clus- 
ter A1689 (z = 0.1832, IStruble & Rood] 1 19991) , a hot, mas- 
sive cluster with a virial mass M vir ~ 1 - 1.5 x 10 15 /r'M w 
i Limousin et al J 1 2007 ; lUmetsu & Broadhursj 120081 : iLemze et al] 



l2009t|Peng et alj200S ) and a n average emission weighted gas tem- 
perature of r ew ~ 10.5 keV {Lemze et alj 120081 : ICavagnolo et alj 
l2009l : iKawaharada et alj[201oF ~ This cluster had previously been 
observed with the CBI1 and thus serves as a useful check on the 
cross-calibration of the two instruments. It also highlights the im- 
proved sensitivity of the CBI2 for SZ observations. To allow direct 
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Figure 6. The distribution of mean rms noise per 4.2 s sample for all the 
CBI2 baselines. The expected value for the nominal system temperature is 
1.9 Jy per sample; the tail to higher values reflects the small fraction of 
antennas with higher than nominal noise. 



comparison of our observations with measurements we also fit the 
combined dataset from the CBI1 and CBI2 to a single isothermal 
beta model. 



5.1 Observations and data 

Observations of A 1689 were carried out with CBI1 and CBI2 over 
the periods 2004 May - June and 2008 January - May respec- 
tively, with the pointing centre at RA(2000) = 13 h ll m 29?5 and 
Dec. (2000) = -01°20'10"0. For both observa tions we adopted a 
simila r observing strategy to that described by lUdomprasert et alj 
( l2004h , whereby any strong correlated ground signal is subtracted 
out using reference fields separated by 8 minutes in right ascension. 
This procedure increases the noise level in the data by a factor of 
i/2 in the case of a single reference field, or \/3/2 in the case 
of averaging over two reference fields. Flagging and calibration 
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Figure 7. Left: The real part of the visibility as a function of baseline for observations of A1689 using the CBI1 (grey errorbars) and CBI2 (black errorbars). 
The short errorbars represent the 1 o~ noise from the variance in the data, and the long errorbars represent the total statistical uncertainty including the intrinsic 
CMB anisotropy. There is about 50 per cent more useable data in the CBI1 observations than in the CBI2 observations. Middle and Right: CLEANed full- 
resolution maps of A1689 using the CBI1 and CBI2 arrays respectively. The rms noises on each map are 5.9 mJy beam -1 for CBI1 and 6.1 mJybeam -1 for 
CBI2. The contours are multiples of 6 mJy beam - 1 . The full-width at half-maximum of the synthesized beam is shown in the bottom left-hand corner of each 
map (8.6 x 8.7 arcmin for CBI1 and 5.3 x 6.0 arcmin for CBI2). The full-width at half-maximum of the primary beam, at v = 31 GHz, is shown as a dark 
circle centered on each map (45. 1 arcmin for CBI1 and 28.2 arcmin for CBI2). 



of the visibility data were performed using CBICAL, a special- 
ist data reduction package designed for use on CBI data. A total 
of 20,259 good visibilities were collected by CBI1 and 13,295 by 
CBI2, where a visibility represents a single 8-minute scan on each 
of the main and trail fields for each of the 78 baselines and 10 fre- 
quency channels. This corresponds to a total equivalent observing 
time (main plus trail fields) of 6.9 hours by CBI1 and 4.5 hours 
by CBI2. The amplitude and phase were calibrated to nightly ob- 
servations of Jupiter and Tau A, or a suitable unresolved source 
when neither primary source was available. The calibration is ul- 
timately tied to the measure d brightness tem perature of Jupiter at 
33 GHz, Tj = 146.6±0.75 K faill et al.l2009l) . and this introduces a 
0.5 per cent calibration uncertainty in the data. In addition to the ab- 
solute flux calibration, short observations of secondary calibrators 
(including 3C 273, 3C 274, 3C 279, J1924-2292 and J2253+1610) 
were used to characterize any possible residual pointing error in 
the experiment. These observations show that the data are consis- 
tent with a residual pointing error of 0.5 arcmin. 

Full-resolution maps of both observations are shown in 
Figure__ and w ere deconvolved using the CLEAN algorithm 
( lHogbom|[T974T) implemented by the APCLN task in AIPS_]. The 
decrement in brightness due to the thermal SZ effect can clearly 
be seen at the centres of both maps. The calibrated visibilities 
were gridded into regularly spaced estimators using the MPI- 
pr ogram , which implements the technique developed by 
iMvers et al.l d2003l) . This significantly reduces the amount of data 
and the size of the covariance matrices that need to be processed 
during the model fitting. 

5.2 Additional sources of error 

5.2.7 Intrinsic CMB anisotropy 

In addition to the instrumental thermal noise and calibration, a sig- 
nificant source of uncertainty for CBI1 SZ observations is the pres- 



ence of the intrinsic CMB anisotropy dUdomprasert et all I2004T) . 
The CMB contributions in the gridded visibilities are correlated 
and therefore must be treated by the construction of a covariance 
matrix. The intrinsic CMB covariances at the positions of the grid- 
ded visibility data are constructed based on a model estimate of 
the CMB power spectrum. MPIGRIDR is used to construct the 
covariance m atrix, with an input power s pectrum generated using 
CMBFASlFl ( Seliak & Zaldarriagall 19961) . assuming a flat ACDM 
cosmology with £1m = 0.3, £1^ = 0.7 and /; = 0.7. The uncertainty 
due to the intrinsic CMB dominates on the largest scales and so has 
the greatest effect on data from the CBI1 array, which has shorter 
minimum baselines than the CBI2 array. Figures[8] and [9] show the 
binned real visibility data as a function of baseline and include the 
estimated uncertainty due to the intrinsic CMB anisotropy. 



5.2.2 Point sources 

Bright point sources can contaminate the SZ signal, appearing as 
positive sources if in the main field or negative if in a reference 
field. Sources that are close to the field centres can have a sig- 
nificant effect on the SZ decrement, especially since there is rel- 
atively low attenuation from the primary beam. Spatial filtering of 
the CBI1 and CBI2 data to angular scales smaller than 10 arcmin 
(corresponding to baselines longer than 300 wavelengths) reveals 
no significant point source fl ux density a b ove th e noise level. Two 
sources were identified by iReese et al.l d2002l> at 30 GHz from 
BIMA and OVRO observations of A 1689. These sources are at po- 
sitions (RA(2000),Dec. (2000)) = (13 h ll m 31?6, -01°19'33f'0) 
and (13 h ll m 30?l,-01°20'37"0), and have integrated flux densi- 
ties of 1.33±0.10mJy and 0.45±0.09mJy respectively. The flux 
density contribution is unlikely to cause significant contamination 
to the CBI data; however, we do subtract these sources from the 
gridded visibility data. Any further residual point source error is 
accounted for in the systematic error estimate for model fitting. 



6 http://www.aips.nrao.edu/ 

7 http://www.cita.utoronto.ca/~myers/ 



http://lambda.gsfc.nasa.gov/toolbox/tb_cmbfast_ov.cfm 
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Table 1. Posterior estimates of the model parameter A7"o and derived parameters F 2 500> ^7500 > ^200 an d ^200 ' or A1689 from fitting to CBI1 and CBI2 data. 
The other model parameters are j8, which is given a strong Gaussian prior of j3 = 0.688 ± 0.013, and r colc , which is given either a strong Gaussian prior of 
^core = 68.4 ±2.1, or a flat uniform prior. Y2S00 and F200 are calculated by integrating the derived Comptonization parameter within projected radii of ^2500 
(200arcsec) and f-200 (600 arc-see) respectively using the strong prior on r corc . I2500 an£ ! ^200 are derived from model fitting to the data with a uniform prior o n 



. A7q is derived from the strong prior on r colc only. Error intervals represent 68 per cent confidence. Priors on r colc and 



to the data with a uniform prior 011 
/3 are from lLaRoque et alj 420061) . 



Parameter 



CBI1 



CBI2 



Combined 



ATo (mK) 


-1.09^ 


-1 23+°- 24 
1 - z -'-0.26 


-1 ?n+ - 18 

1 - zu -0.18 


F 25 oo (10- 10 sr) 


1.79l»i 


2 02+ 42 


1 95+ 033 

1 '"-0.28 


y 2 %(io- 10 sr) 


1 1 7 +0.56 
• -Q.45 




ioo -0.38 


K 2Q0 (10- 10 sr) 


7.00+}^ 


8.12^ 


7 71+1.48 
'•' 1 -1.16 


T|g(10- 10 sr) 


3.65+?f 


5.79^ 


s 7S+1.68 
■'•■"-1.31 



5.2.3 The Kinematic SZ effect 

The non-zero peculiar velocities of galaxy clusters cause a sec- 
ondary distortion in the C MB frequency spectr um known as the 
kinematic SZ (KSZ) effect. iBenson et al.l l l2003h estimated the pe- 
culiar velocity of A1689 to be v pec = +170+^ ±750kms _1 using 
the SuZIE II experiment (where the first quoted errors are statisti- 
cal and the second systematics). Clearly the uncertainties in this 
measurement dominate, however we can estimate the uncertainty 
introduced by the KSZ effect by assuming a typical line-of-sight 
pecul i ar velocity of ±300kms~ 1 1 Watkins 1 19971 ; iGiovanelli et alj 



Il998t iDale et al.ll 19991 : IColberg et al .11 200(3) . At an observing fre 
quency of 3 1 GHz the non-relativistic KSZ effect signal for a 
10.5 keV cluster is ±2.6 per cent of the thermal SZ signal. Further 
relativistic corrections to the KS Z signal have been calculated by 
iNozawa. Itoh & Kohvam j < fl998h and change the error by only 0.2- 
0.3 percent of the thermal SZ signal for high temperature clusters. 



5.3 Model fitting 

We fit an analytical model to the visibility data in order to derive 
properties of the cluster that can easily be compared with the liter- 
ature values. We assume that the observed thermal SZ effect has 
circular symmetry and we model the change in brightness tem- 
perature using the single isothermal jS model given by equation 
[2] Model fitting is done in visibility space where the instrumental 
noise covariance matrix is diagonal (although radio sources and the 
CMB introduce off-diagonal elements to the covariance matrix). 
Imple mentation of the model fitting i s performed using MULTI- 
NEST ( iFeroz. Hobson & Bri dges 2009), a po werful Bayesian opti- 
mizer that uses the Nested Sampling method dSkillindl2004l) . This 
program returns a weighted sampled posterior probability distribu- 
tion for each of the model parameters, which can then be marginal- 
ized over in order to obtain estimates of the derived cluster proper- 
ties. We also introduce a calibration error as a nuisance parameter 
with a Gaussian prior that accounts for a total systematic error of 
5 per cent (one sigma), and which is later marginalized over. 



5.4 Results 

TableQ] shows the priors and estimates of model parameters from 
fitting separately to CBI1 and CBI2 data, and then jointly to both. 
The (11, v) coverage of the CBI arrays means that r cole and j8 are not 
individually well constrained by the data and so we apply Gaus- 
sian priors to these p arameters based on the values measured by 
lLaRoque et alj | |2006|) from a combined X-ray and SZ analysis. 
The estimated values of the central SZ decrement for the CBI1 and 



CBI2 arrays are consistent within the errors. Although the error 
bars on the parameters (which are derived from the posterior prob- 
ability distributions) are similar for CBI1 and CBI2, the joint fit is 
dominated by the CBI2 data. This is because the CBI1 data have 
significant off-diagonal terms in the covariance matrix due to the 
intrinsic CMB fluctuations, which do not integrate down with the 
addition of more data. Our joint fitting to both data sets fully takes 
this effect into account, but it means one cannot simply take the 
weighted mean of the individual CBI1 and CBI2 parameter esti- 
mates. 

From the model we can calculate the total integrated y- 
parameter for a given projected aperture. This quantity is a measure 
of the cluster's total thermal energy contained within the radius of 
integration, and is given by 



Y = 



I ydfl = 2k Jyrdr, 



(4) 



where r is the projected angular radius from the centroid. Y is the 
most interesting observable parameter in terms of relating mea- 
sured SZ signals to the intrinsic cluster properties such as mass. 
In order to compare with values given in the literature we calculate 
^2500 within a radius of 200 a rcsec , which correspond s roughly to 
an overdensity radius of r^nn. iBonamente et alj i2008l) used obser- 
vations of A 1689 with the BIMA and OVRO arrays to obtain an es- 



timate of Y' 



2500 



1. 



f0.49 
^-0.38 



x 10 1 sr with (-2500 = 1 96 ± 8 arcsec , 



and iLiao et all fcoiOh used AMIBA observations to obtain an es- 



timate of *2500 = 3. l + [ 3 x 10~ lu sr with r 2 soo = 215i[9 arcsec. 
Figure[8] shows the posterior probability distribution for estimates 
of Y from CBI1 and CBI2 observations. Our estimate of ^SOO = 
1.95^Q 2g x 10~ 10 sr is consistent with these results, within the er- 
rors. 

We also calculate Y out to a larger projected radius of 
600 arcsec, wh i ch is approximately the value of r200 quoted by 
lLaRoque et al.l I2OO6), in order to demonstrate that the CBI1 and 
CBI2 arrays are capable of measuring the thermal SZ effect inte- 
grated out to large physical radii. The value of Y at 7-200 was cal- 
culated using the strong prior on r core 

used above, giving Y200 = 
7.71^j jg x 10~ 10 sr. This significantly larger value than J^goo re- 
flects the fact that a large fraction of the integrated cluster pressure 
is outside the core region probed by higher resolution experiments. 
The fit to the CBI visibility data and the posterior probability dis- 
tributions for 1200 and ^2500 are shown in Figure [8] 

We have used prior values of j8 and r core derived from X-ray 
measurements to estimate A7o, as our SZ data by themselves do 
not put strong constraints on these parameters. However, it is not 
necessary to accurately constrain all the parameters of the /3 model 
in order to make a good measurement of Y. We therefore also cal- 
culate Y using a flat prior on r core , i.e., assuming we have no 



16 , 
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Figure 8. Results from the fit of an isothermal beta model to the combined CBI1 and CBI2 data, using the priors on j0 and fcore 

from lLaRoque et alj feOOd) 

(J3 =0.688 ±.013, r coie — 48.4 ±2.1 arcsec). Left: The real part of the visibility for the combined CBI data set as a function of baseline distance. The short 
errorbars represent the 1 o~ noise from the variance in the data, and the long errorbars represent the total statistical uncertainty including the intrinsic CMB 
anisotropy. The grey scale represents the 68.3, 95.4 and 99.7 per cent confidence intervals of the fitted model. Middle and Right: The estimated posterior 
probability distributions for the integrated Comptonization parameters i2500 an d ^200- 
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Figure 9. (Top) As Figure [8]but with a uniform prior on r corc , showing that reasonable constraints can be placed on Y even without detailed prior knowledge 
of the cluster core size. (Bottom) The joint constraint on A7q and r corc , showing that they are individually poorly constrained. 



knowledge of this parameter, while maintaining the strong prior 
on j6, which typically does not vary significantly from a value of 
j6 ~ 2/3 between different clusters. Figure [9] shows the fit to the 
CBI visibility data, the probability distributions for Y^qq and Y^L., 
and the likelihood contours for the free parameters r cole and ATq. 
Although r core and A7y are individually very poorly constrained, 
the fits for K flat are only slightly worse than when r core is strongly 
constrained. The estimates of F flat are also significantly improved 
when fitting to the combined CBI1 and CBI2 data sets. This is to 
be expected since Y is proportional to the SZ total flux density (i.e., 
the zero-spacing visibility), which is measured almost directly by 
the CBI short baselines. 



commercial off-the-shelf antennas with custom-made secondary 
mirrors. The fabrication techniques were designed to be cheap and 
easy to reproduce and could be used as a cost-effective method for 
producing a large number of antennas in this size range. The new 
antennas have been demonstrated to perform as specified, and us- 
ing SZ observations of the cluster A 1689 we have shown that the 
inter-calibration of CBI2 and CBI1 is good, that the upgraded array 
met its design sensitivity, and that CBI2 and combined CBI1 plus 
CBI2 data can be used to constrain the integrated SZ Comptoniza- 
tion parameter out to large radii 0'2oo) both with and without prior 
information on the cluster size. 



6 CONCLUSIONS 

We have described an upgrade to the Cosmic Background Imager in 
which the original thirteen 0.9 m antennas were replaced with new 
1.4 m antennas. The upgrade was achieved by using inexpensive, 
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